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ABSTRACT 

The static electrical conductivity of non-ideal, dense, partially ionized he- 
lium plasma was calculated over a wide range of plasma parameters: tempera- 
tures 1 • 10 4 K < T < 1 ■ 10 5 K and mass density 1 x l(T 6 g/cm 3 < p < 2g/cm 3 . 
Calculations of electrical conductivity of plasma for the considered range of 
plasma parameters arc of interest for DB white dwarf atmospheres with ef- 
fective temperatures 1 • 10 4 K < T eff < 3 • 10 4 K. 

Electrical conductivity of plasma was calculated by using the modified 
random phase approximation and semiclassical method, adapted for the case 
of dense, partially ionized plasma. The results were compared with the unique 
existing experimental data, including the results related to the region of dense 
plasmas. In spite of low accuracy of the experimental data, the existing agree- 
ment with them indicates that results obtained in this paper are correct. 
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1 INTRODUCTION 



DB white dwarf atmospheres belong to the class of astrophysical objects wh i ch have been 



investigated for a long time and from various aspects (IBues 



1970 



Koester 



1980 



Standi et al. 
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19931 ). In the previous period, the contribution of the authors of this paper was related to 



Mihailov et al 



Ignjatovic et al. 



(2003) - to the 



the research of optical properties of DB white dwarf atmospheres wi thin the range of average 
effecti ve temp eratures 1 ■ 10 4 K < T P ft < 2 ■ 10 4 K. So, the papers of |]VIihajlov &: Dimitrijevi£ 
( 119921 ) and of iMihajlov et al.l (11994 Il995[ ) were dedicated to the research_ of con tinual ab 
sorption in the optical part of EM spectra, the paper of 
research of chemi-ionization/recombination processes and of 
investigation of continual absorption in the VUV region of EM spectra. 

Recently, the transport properties of helium plasmas, characteristic of some DB white 
dwarf atmospheres, attracted the authors' attention, first of all the electrical conductivity. 
Namely, the data on electrical conductivity of plasma of stars with a magnetic field or mov- 



(120091 ) - to the 



i ng in t he magnetic fie l d of t he ot her component in a 



((20091), 



Potter fc Tout 



2009) and 



j mary 



Rodriguez-Gil et al. 



system (see e.g. 



Zhang et al. 



(120091 )) could be of significant in- 



terest, since they are useful for the study of thermal evolution of such objects (cooling, 
nucl ear burni n g of a ccret ed matter) and the i nvesti gation of their magnetic fields. For exam- 
ple, iKopeckyl ( 119701 ) and iKopecky fc Kotrd ( 119731 ) studied electrical conductivity for stars 
of various spectral types, i n order to investigate the magnetohydrodynamic differences in 



their atmospheres. Recently, 



Mazevet et al. 



(120071 ) investigated He conductivity in cool white 
dwarf atmospheres, since the possibility of using these stars fo r dating stellar popu lations 



has generated a renewed interest in modelling their cooling rate ( iFontaine et al. 



200lh. Also. 



Baiko Sz Yakovlev 



the tra nsport processes occurring in the cores of white dwarfs (see e.g. 
( 1 19951 ) and numerous references therein) have been considered. Moreover, electrical conduc- 
tivity was particularly investigated for solar plasma, since it is of interest for consideration of 
various processes in the observed atmospheric layers, like the relation between magnetic field 
and convection, th e question of m agnetic field dissipation and the ener gy released by s uch 



processes (see e. g. 



Kopeckyl ( 119701 ) and references therein) . For example iFeldmanl ( 119931 ) in- 



vestigated the role of electrical conductivity in the construc t ion o 



upper Solar atmosphere, and 



Kazeminezhad &: Goodman! ( 120061 ) considered the electrical 



a theoretical model of the 



conductivity of solar plasma for magnetohydrodynamic simulations of the solar chromo- 
spheric dynamo. Given that electrical conductivity plays an analogous role in other stars as 
well, it is of interest to investigate its significance, to adapt the methods for research into 
stellar plasma conditions and to provide the needed data. 

An additional interest for data on electrical conductivity in white dwarf a tmosp heres 
may be stimulated by the search for extra-solar planets. Namely iJianke et al.l (119981 ) have 
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shown that a planetary core in orbit around a white dwarf may reveal its presence through 
its interaction with the magnetosphere of the white dwarf. Such an interaction will generate 



Jianke et al. 



electri cal currents that will directly heat the atmosphere near its magnetic poles. 
(119981 ) emphasize that this heating may be detected within the optical wavelength range as 
H a emission. For investigation and modelling of mentioned electrical currents, the data on 
electrical conductivity in white dwarf atmospheres will be useful. 



One of the most frequently used approximations for considerat ion of transport properties 



ofdij 


Terent plasmas is t 


he approximation of " fully ionized plasma" ( 


Spitzei 


1962 


Radke et al. 


1976 




Adamvan et al. 


1980 


• 


Kurilenkov 


1984 




Ropke & Redmer 


1989 




Diuric et al. 


1991: 


Nurekenov et al. 


1997; 


Zaika et al. 


2000 


Esser et al. 


2003 


). It was shown that the electri- 



cal conductivity of fully ionized plas mas can be success 



random-phase approximation (RPA) (iDjuric et al. 



1991 



ully calculated using t he modified 



Adamyan et al 



1994a 



bf) in the re- 



gion of strong and moderate non-ideality, while the weakl y non-ideal plasmas were success 



fully treated within the semiclassical approximation (SC) (jMihajlov et al 



1993 



Viteletal 



200ll ). In practice, even the plasmas with a significan t neutral component are treated as fully 



ioniz e d in order to simp l ify the considered problems, ( iRopke fc Redmerlll989i : lEsser fc Ropke 



1998 



Zaika et al 



2000 



Esser et al. 



2003). However, our preliminary estimates have shown 



that such an approac h is not applic able for the helium plasmas of DB white dwarf atmo- 



spheres described in (IKoester 
neglected. 



19801 ). where the influence of neutral component cannot be 



Therefore, an adequate method for calculations of electrical conductivity of dense, par- 



tially ionized helium plasmas is de veloped in t 



alization of methods developed in 



Djuric et al 



l is pap er. T 
(I19911 ) and 



lis method rep r esents a gener 



Mihajlov et al. 



(119931 ) . namely 



modified RPA and SC methods, and gives a possibility to estimate the real contribution of 
the neutral component to the static electrical conductivity of the considered helium plasmas 
within a wide range of mass densities (p) and temperatures (T). 



The calculations were performed for helium plasma in the state of local thermodynamical 
equilibrium with given p and T for 1 ■ 10 4 K < T < 1- 10 5 K and 1 x 10 _6 g/cm 3 < p < 2g/cm 3 



he o b tained results are compared with the corresp onding experimental data (IMintsev et al. 



(I1980h : 



Ternovoi et al 



(|2002|); 



Shilkin et al 



(120031 ) ). Fo r the calculatio ns of plasma charac 



teristics of DB white dwarf atmospheres the data from (lKoesterlll980l ) were used. 
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2 THEORY 

2.1 The plasma electrical conductivity 

On the basis of the previous paper of 



Adamyan et al. 



f|l980h . a modified RPA method 



for ca lcula tion of static conductivi ty of fully ionized plasma was developed in iDjuric et al. 



(1199lh and 



Adamyan et al. 



(1l994af l. The principal role in this method belongs to the formula 
for energy- dependent electron-electron (ee) and electron- ion (ei) relaxation times t ee;e i(E) = 
tfveii where E is the energy of single electron state, determined as a sum over the Matsubara 
frequencies by using the methods of Green function theory. This method is especially suitable 
for calculation of electrical conductivity of dense non-ideal plasmas with electron density (N e ) 
larger than 10 17 cm -3 . In the region N e < 10 17 cm -3 the static cond uctivity of fully i onize d 
plasmas can be determined well using the SC method developed in iMihajlov et al.l (119931 ). 
which is also based on electron relaxation time t ee . ei (E) = tf e c ei . It is important that SC 
method gives practically the same results as RPA method in a wider region of the electron 
densities around the value of N e = 10 17 cm~ 3 . SC method was tested from this aspect in 



Vitel et al. 



from 



(2001), wher e it was experiment ally verified through comparison with the results 



Spitzerl (119621 ) and iKurilenkovl (1 19841 ). just for the helium plasmas. 



However, as it was already mentioned, the helium plasma of the considere d DB white 



leres contains significant neutral atom component, as follows from lBuesl (1l970h 



dwa rf atmosp 

and iKoesteii (Il980l ). Because of that, we will start here from the fact that in both RPA and 

SC methods effective electron relaxation times t ee - e i(E) can be expressed as 

1 



(1) 



where E is the electron energy, and u ee . ei (E) - the corresponding total electron-electron and 

electron-ion collision frequency. This gives a possibility to generalize the modified RPA and 

SC methods for the case of partially ionized plasmas, replacing l/t ee . e i(E) by u e] t t(E) 

1 



u e;tot (E) = v ee . ei {E) + u ea (E) 



Vea(E), 



(2) 



(E) 

where v ea [E) is an effective electron-atom collision frequency. Consequently, the basic RPA 
and SC expressions for the static electrical conductivity o"o transform to the corresponding 
Frost- like expressions 

1 df FD (E) 



4e f°° 
*°=3Wo E - W{E) - 



+ Vea{E) 



dE 



-dE, 



(3) 



(E) 

where m and e are the mass and the modulus of charge of the electron, w(E) is the density 
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of the single electron states in the energy space, fFo(E) = fpD(E; T, N e ) is the Fermi-Dirac 
distribution function for given N e and temperature T, and t ee;ei (E) = t^ e f or t ee]ei (E) = 



On the basis o: 



atmospheres (iBues 



composition and te mperature data of the considered DB white dwarf 



1970 



Koester 



1980) , one can note that only four components of these 



atmospheres are important for determination of a^. free electrons, He + and H e 2+ ions, and 
helium atoms. In accordance with this, the expression for 1/t^f from 
can be presented in the form 



Djuric et al. 



(1991) 



1 _ ATimN e e 4 k B T 



VSmE/h ^ 




E 

j=l,2 



NijeKq) 



(4) 



where h and k B are the Planck and Boltzmann constants respectively, Z e ■ e and • e - 
the charges of electron and He + (ls) and He 2+ ions (Z e = —1, Z a = 1, Z i2 = 2), Nn and 
- the corresponding helium ion densities, n e;n (g), Ha. n (q) and n i2;n (9) - the electron 
and ion polarization operators, e n (q) - the dielectric function, n = 0, ±1, ±2..., and summa- 
tion is extended over all the Matsubara frequencies Q n = 2irnkBT /h. The detai l ed ex pres- 



sions for polarization op erators and dielectric functions are given in 



Adamvan et al. 



Djuric et al. 



f|199lh and 



(119941 . 

The SC expression for <To is applied to the outer layers of the considered DB white dwarf 
atmospheres where A^ e < 10 17 cm -3 and the presence of H e 2+ ions can be n eglected. Because 



of that the expression for l/tf e c ei , in accordance with 



Mihajlov et al.l (119931 ). can be taken in 



the form 



1 (2m) 1 /2 jB 3/2 



Xee 2e*Z a N e ln(l + A 2 )V2 



A; 



IE 

Zne 2 



(5) 



where k is Boltzmann constant and Zn = 1. The correction factor 1/Xee — ^/Xee{Zn,T) is 
determined within SC method, while r ci ~ (A7re 2 N e /kT)~ 1 ^ 2 is the corresponding screening 
length which is an external param eter of the theory. Th e values of 1/Xee f° r %n = 1 an d 



Mihajlov et al.l (I1993L Here, value r 



is taken as 



10 4 K ^ T ^ lO 5 ^ are taken from 
th e screening leng t h, wh ere the ion neutrality radius r n;i is given by the expressions obtained 



m 



Mihajlov et al. 



(120091 ) . Such a choice of the screening length provides overlapping SC and 
RPA values of a in a wider region of the electron densities around the value of N e = 10 17 



cm 



Finally, the electron- atom collision frequency v ea (E) in Eq. ([3]) is given here by the known 
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expression 

u ea (E) = N a -v(E)-Q t ; a (E), (6) 

where N a is the He(ls 2 ) atom density, v(E) = (2E /m) 1 / 2 - relative electro n-atom velocity, 
and Q l L (E)- the transport cross-section for the elastic e— He(ls 2 ) scattering (IMott fc Massey 
1970). 



2.2 The electron-atom transport cross-section 



The exact quantum-mechanical calculation of the transport cross-section for elastic electron- 
atom collisions is a very hard problem in itself. However, in the case of e — He(ls 2 ) scattering 
the problem becomes easier within the considered temperature range where all non-elastic 
collision processes can be neglected. Namely, in this case it can be treated as scattering of 
electrons in an adequately chosen model potential U(r), where r is the distance between the 
electron and the nucleus of He(ls 2 ) atom. Thus, the electron movement is described by a 



wave function \I/(r, 9, 



Xi\r) 



Y, 



lm \ 



where Y\ r 



(f)) is a spherical harmonic function 
r 

of degree I and order m, while the function xi{ r ) satisfies the radial Schrodinger equation 

ld 2 r „, , 1(1 + 1 



2 + U(r) + 



X(r) = Ex{r), 



(7) 



2 dr 2 ' ~ ' ' 2r 2 

given in atomic units. Following the previous paper of llgnjatovic fc Mihajlovi (jl997l ). we take 
the model potential U(r) in the form 

Z q Z — q 



U{r) 



U (r) - 

r r + r 

U m (r) = ar 2 + br + c, 
U as (r) 



< r < r,- 



a 



r 

r,i < r < r f 
r j < r < oo. 



(8) 



2(r 2 + h 2 Y 

where Z = 2 is the charge of the nucleus of a helium atom, q and Z — q describe the redis- 
tribution of electrons in the Is shell of the helium atom (0 < q < 2), a is the polarizability 
of He(ls 2 ) atom, and ro has the meaning of the average radius of atom. Parameters a, b, 
c and h are determined from the conditions of continuity and smoothness of the potential 
U(r) at the points and 77 where < r, < ro and rj > r^. The values of the mentioned 
quantities are given in Table (IIJ). 



The equation (J7|) is solved by the partial wave method described in llgnjatovic fc Mihajlov 



(Il997l ). One obtains as the result the phase shifts 81(E) for I = 0, 1,2, where each 81(E) 



corresponds to the partial wave with a given orbital number I. This way of solving Eq. (J7J) 
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Table 1. The parameter values of the model potential U(r) from Eq. ||5} in atomic units 



z 


'J 


a b 


c ah 


ri rj r 


2 


0.7 


-0.59597652 2.2545472 


-2.19405731 1.384 0.01 


0.73 1.75 0.9 



has advantages since, apart from the transport cross-section 

Q tr (E) = JL- YJLI + l)stn 2 (5 l - 5 W ), (9) 

1=0 

it allows determination of the elastic cross-section, as well as cross-sections of higher order 
(viscosity cross-section, etc.). As one can see in Figure HJ there is an excellent agreement 
between the values of Qt r {E) calculated using the potential (jSJ) and various experimental 
and theoretical data of other authors. 



3 RESULTS AND DISCUSSION 

In order to apply our results to the study of DB white dwarf atmosphere plasma properties, 
helium plasmas with electron (N e ) and atom (N a ) densi ties and temp eratures (T), charac- 
teristic of atmosphere models presented in the literature (]Koesterlll980l ). are considered here. 
So, the behaviour of p and T for models with the logarithm of surface gravity log g =8 and 
effective temperature T e ff =12000K, 20000K and 30000K is shown in Fig. [2] as a function 
of Rosseland opacity r. As one can see, these atmospheres contain layers of dense helium 
plasma. In order to cover the considered plasma parameter range reliably, we tested our 
method for calculation of the plasma electrical conductivity within a wider range of mass 
density 1 x lCr 6 g/cm 3 < p < 2g/cm 3 and temperature 1 • 10 4 K < T < 1 • 10 5 K. 

The influence of neutral atoms on the electrical conductivity of helium plasma is shown 
in Fig. El In this figure the electrical conductivities for T = 15000, 20000 and 25000K 
are given as functions of mass density p. The range between the two vertical dashed lines 
corresponds to the conditions in the considered DB white dwarf atmospheres. Two groups 
of curves, calculated using the expression are presented in this figure: a) the dashed 
ones, obtained by neglecting the influence of atoms, i.e. with v ea = 0; b) the full-line curves 
calculated with the influence of atoms included, i.e. with v ea given by Eq. (jHJ). First, one 
should note that the behaviour of these two groups of curves is qualitatively different: the 
first one increases constantly with the increase of p, while the other group of curves decreases, 
reaches a minimum, and then starts to increase with the increase of p. One could explain 
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such behaviour of the electrical conductivity by the pressure ionization. This figure also 
clearly shows when the considered plasma can be treated as "fully ionized". 

In Fig. H] we compare our values of the helium plasma conductivity, shown by full curves 
for T = 15000 K, 20000 K and 25000 K within the region 5 • 10~ 4 g/ cm3 < P < 2g/cm 3 , with 
the existing experimental data. Let us note that these experimental results are uniquely 
available for comparison so that, in spite of their low accuracy, the agreement with them 
gives the only possible indication that our results are correct. 

Within the region p < 0. 65 g/cm 3 , i.e. to the left of t he vertical lin e in F ig. HJ there 



are experimental results from 



Ternovoi et al. 



(11999|) (o) and 



Shilkin et al. 



( 120031 ) (V) where 



the temperature was determined with an error of less than 20%, which are related to the 



the vertical lin e in Fi g. HJ 
Ternovoi et al.l (120021 ) for 



temperature range 20000K-25000K. For p > 0.65 g/cm 3 , i.e. right o 
are shown several values of the plasma conductivity, obtained by 
the temperature range 15000K-25000K. These experimental values are obtained with an 
experimental error ~ 50% and can be treated only as characteristic of this temperature 
region as a whole. These data at least indicate that our results for p > 0.65 g/cm 3 lie in the 
correct domain of the electrical conductivity values. 

The developed method was then applied to calculation of plasma electrical conductivity 
for the models of DB white dwarf atmospheres presented in Fig. [2j The results of the 
calculations are shown in Fig. [51 First, let us note a regular behaviour of the static electrical 
conductivity which one should expect considering the characteristics of DB white dwarf 
atmospheres. Further, the electrical conductivity profiles presented in this figure show that, 
for the considered DB white dwarf models, plasma electrical conductivity changes over the 
domain of values where our results agree with the experimental ones (see Fig. HJ). This 
indicates that the theoretical apparatus presented here may be adequate to be used for 
investigation of DB white dwarfs in the magnetic field of their partners in binary systems 
and magnetic white dwarfs. 

In order to provide possibility for direct applications of our results to different theoretical 
investigations the values of static electrical conductivity <7o of helium plasma in a wide range 
of p and T are given in Table J2J This table cov e rs pla sma conditions for all models of DB 



white dwarf atmospheres presented in 



Koesterl (119800 and the corresponding values of o"o 



were determined for the helium plasmas in the state of local thermodynamical equilibrium 
with given p and T. 

The method developed in this paper represents also a powerful tool for research into 
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Table 2. Static electrical conductivity of helium plasma <tq [l/(Ora)] 



T[K] 










r / 3i 










5.00E-07 


1.00E-06 


5.00E-06 


1.00E-05 


5.00E-05 


1.00E-04 


5.00E-04 


1.00E-03 


5.00E-03 


oUUU 


o.ooizj-ruu 


^.oy-Cj-t-uu 


o 1 ^Tr_i_nn 

Z. 10-Cj-tUU 


i ^np i nn 
l .oun*-|-uu 


u.yvjiij-ui 


A sqt? m 


9 9np ni 




a 09 


yuuu 








l . loej-\-\j ± 


0. 1 ZJZj^UU 


o.DU-cj-ruu 




i in!? i on 


nnr ril 


1UUUU 






7ATT-L.ni 






1 79F-l_ri1 




o. I orj-|"UU 




12000 


9.36E+02 


8.04E+02 


5.09E+02 


4.05E+02 


2.19E+02 


1.66E+02 


8.35E+01 


6.08E+01 


2.92E+01 


14000 


2.27E+03 


2.08E+03 


1.62E+03 


1.38E+03 


8.87E+02 


6.97E+02 


3.99E+02 


3.00E+02 


1.28E+02 


16000 


3.64E+03 


3.55E+03 


3.15E+03 


2.87E+03 


1.45E+03 


1.18E+03 


6.85E+02 


5.28E+02 


2.75E+02 


18000 


4.90E+03 


4.97E+03 


4.76E+03 


3.79E+03 


2.87E+03 


2.46E+03 


1.58E+03 


1.26E+03 


7.07E+02 


20000 


6.08E+03 


6.22E+03 


6.32E+03 


5.61E+03 


4.99E+03 


4.46E+03 


3.21E+03 


2.53E+03 


1.58E+03 


25000 






9.95E+03 


1.02E+04 


9.90E+03 


9.90E+03 


8.62E+03 


7.74E+03 


5.95E+03 


30000 






1.33E+04 


1.40E+04 


1.56E+04 


1.60E+04 


1.57E+04 


1.50E+04 


1.21E+04 


35000 






1.65E+04 


1.74E+04 


2.02E+04 


2.15E+04 


2.34E+04 


2.35E+04 


2.15E+04 


40000 






1.97E+04 


2.09E+04 


2.43E+04 


2.61E+04 


3.06E+04 


3.18E+04 


3.22E+04 


45000 






2.31E+04 


2.45E+04 


2.84E+04 


3.06E+04 


3.69E+04 


3.93E+04 


4.36E+04 


55000 






3.04E+04 


3.23E+04 


3.70E+04 


3.98E+04 


4.81E+04 


5.26E+04 


6.40E+04 


65000 






3.81E+04 


4.08E+04 


4.65E+04 


4.98E+04 


6.00E+04 


6.56E+04 


8.20E+04 


75000 






4.58E+04 


4.93E+04 


5.66E+04 


6.06E+04 


7.26E+04 


7.92E+04 


9.84E+04 



white dwarfs with different atmospheric compositions (DA, DC etc.), and for investigation 
of some other stars (M-type red dwarfs, Sun etc.). Finally, the presented method provides 
a basis for the development of methods to describe other transport characteristics which 
are important for the study of all mentioned astrophysical objects, such as the electronic 
thermo-conductivity in the star atmosphere layers with large electron density, electrical 
conductivity in the presence of strong magnetic fields and dynamic (high frequency) electrical 
conductivity. 
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Figure 1. Momentum trans fer cros s secti on Q tr as a functio n of energy E. C alculated data (— line ) for helium tog ether 
with results of o t her a uthor s iNesbetl il979l) (- - das hed linel.lFon et al] dl98lij (• • • dotted linel , iReeister et al.l dl980h (A), 
ICrompton et alj dl970T) (*), iBrunger et al.l dl992f) (o1. iMillov fc Cromptonl <1977l) (□). 
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Figure 2. DB white dwarf atmosphere models with log 3 = 8 and T e j^ = 12000K (full curve), T e yy=20000K (dashed curve) 
and T e// =30000K (dotted curve) from lKoesteil dl98(J) : (a) The mass densities; (b) The temperatures, as functions of Rosseland 
opacity r. 
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Figure 3. Static electrical conductivity <rrj of dense He plasmas as a function of mass density p (full curves), compared to the 
Coulomb part of conductivity (dashed curves). The area left to the vertical dashed line marks the region which is of interest 
for DB white dwarfs. 
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Figure 4. Static electrical conductivity no of helium plasma for various temperatures as a function of the mass density p. Full 
line - the calculations based on the expressions (0-(|6]l; right of the vertical dotted line there is a region of extremely dense 
plasmas where one s hould treat the pres ented calculation as an extrapolation, v - IShilkin et al.l l|2003l ) 20000-23000K; ► - 
iTernovoi et~aH i2002h 15000-25000K; Q - iMintsev et~aH il980h 20000-25000K 
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Figure 5. Static electrical conductivity <ro as a function of the logarithm of Rosscland opacity t for DB white dwarf atmosphere 
models with log g = 8 and T e// =12000K (full curve), T e// =20000K (dashed curve) and T e// =30000K (dotted curve) 



